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^j _ in systems which obey generalized statistics, namely fractional exclusion statistics and 

Gentile's statistics. The transition rates are derived in the framework of canonical 
ensembles. This approach introduces a tool for describing interacting fermionic and 
bosonic systems in non-equilibrium as ideal FES systems, in a computationally efficient 
manner. The two types of statistics are analyzed comparatively, indicating their 

—i- \ intrinsic thermodynamic differences and revealing key aspects related to the species 

"sj" ' size. 
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1. Introduction 

Two of the best known generalizations of the Bose and Fermi statistics are the Gentile's 
statistics (GS) jH |2j [3] and the fractional exclusion statistics (FES) [3]. A method for 
the stochastic simulation of the time evolution of Bose and Fermi gases was proposed 
by Guastella et al. in Ref. j5]. In this paper we present a method for the simulation of 
the more general GS and FES systems, which have as limit cases the Fermi and Bose 
statistics. As we shall see below, the GS and the FES are apparently closely related. 
Nevertheless they give totally different thermodynamic results. 

The GS describes systems of non-interacting particles in which each single particle 
state can be occupied by a maximum of £ particles; the case £ = 1 corresponds to the 
Fermi statistics and the case £ = oo corresponds to the Bose statistics. 

The definition of FES is somewhat more complicated. A FES system is formed in 
general of several subsystems, called species. Each species consists of a finite number 
of particles in a finite dimensional space, spanned by a set of single-particle quantum 
numbers. They are defined by coarse-graining the set of single-particle quantum numbers 
for each type of particles in the system. Each species has a certain number of available 
single-particle states which may depend, in principle, on the number of particles in any 
of the species of the system. This dependence is known as the "statistical interaction" 
between the FES particles and limits the number of particles that can be accommodated 
in each species, similarly as in the GS. 

Let us number the species by i = 0, 1, . . . and denote by Ni and Gi the number 
of particles and the number of available single-particle states of species i, respectively. 
The statistical interaction is expressed by the change in the number of available single- 
particle states at the change of the particle numbers. For example a change SN{ of 
Ni will change Gj by —aji5Ni, for any j. The FES parameters, CKy, depend on the 
way the system is divided into species - different divisions lead to different a's [6]. For 
non-interacting bosons and fermions, a^ = and a^ = 8ij, respectively, for any i and 

3- 

Now the connection between GS and FES is obvious. If we would have two similar 
types of systems, one described by GS, of parameter £, and the other one described 
by FES, of parameters aiy = SijOt = (%/£, we may coarse-grain the sets of quantum 
numbers that describe these systems and obtain the "species" i — 0, 1, . . ., each with Ni 
particles and Gi available states. If we denote by G\ the number of available states 
when Ni = 0, for any i, then the maximum number of particles that we can put into 
the species i in both gases is the same, namely t;G\ = G\ /a. Based on this single 
observation we could say that GS is pure FES with species formed of single-particle 
states. 

Nevertheless, the thermodynamics of FES and GS systems are different. In general, 
FES is applied to species of large numbers of particles and states and cannot be applied 
to single states, like GS. If one would attempt to make the FES species smaller and 
smaller, not only that the thermodynamic results deviate from those obtained with 
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large species, as we shall show below, but the very definition of FES becomes unclear. 
On the other hand, if we apply both, FES and GS to particle species-say we have 
the species {N i: Cri}i=o,i,.. -the number of microconfigurations in which the particles can 
arrange themselves is not the same in FES [U [7J El |9j [10] and GS [U |2j [3] systems. 
Therefore the thermodynamic results must be different. 

The difference between the FES and the GS is very clearly evidenced in systems of 
constant density of states (DOS). It is well known that Bose, Fermi and, in general, FES 
systems of constant DOS are thermodynamically equivalent, i.e. the specific heat and 
entropy are independent of the statistics [HJ [121 031 CHI US [16] . The microscopic reason 
of this equivalence lies in the fact that one can realize a one-to-one mapping between 
microconfigurations of particles of the same excitation energies across all statistics - 
by the "excitation energy" of a microconfiguration we mean the difference between the 
energy of the microconfiguration and the lowest energy of the system [13] . As we shall 
see in the following, the GS does not fall into this equivalence class. For a GS of constant 
DOS, the specific heat, and therefore also the entropy, depend on the parameter £. 

As stated above, we propose in this paper a method to stochastically simulate the 
time evolution of GS and FES systems. The paper is organized as follows. For the clarity 
of the exposition, we start in Section [2] by briefly reviewing the method of Guastella 
et al.. Then we extend the method to the FES and GS gases in Section 3. Finally, in 
Section 4 the derived transition probabilities associated with the Markov chain that we 
introduce are used in simulations of FES and GS systems, the results are tested against 
corresponding analytical data and particularities of the two statistics are commented 
upon. 

2. Simulation of Bose and Fermi ideal gases 

Let us consider a system of ideal bosons or fermions in contact with a heat bath, at 
temperature T. The single particle states of the system, denoted by \i), have the energies 
6j (i — 0, 1, . . .); the energy levels may be degenerated. 

The dynamics of the system is described in terms of a time dependent Markov chain 
in the microconfiguration space. Each microconfiguration is a set {^i}i=o,i,... 5 where rii is 
the particle population of the state \i). For fluency we shall call the microconfigurations 
simply states, which must not be confused with the single-particle states. The system 
evolves in time from one state to another neighboring state by random jumps of particles 
between the single-particle states. The transition rates between two states will be 
denoted by T{ n >.y t { ni y, where the notation convention is that the first subscript refers 
to the final state and the second subscript refers to the initial state. Assuming that in 
the dynamics of the system the particles jump between single particle states one at a 
time, we shall consider that two states are neighbors if they differ by the exchange of 
only one particle between two single particle states. Therefore, we shall use the simpler 
notation Tji(rij,ni) to refer to the transition rate of a particle from the state \i), of n* 
particles, to the state \j), of rij particles. 
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If the particles are bosons, the Fermi golden rule imposes in general that Tji(nj,rii) 
is proportional to (rij + l)rij. For fermions, the Pauli exclusion principle gives a transition 
rate Tji oc (1 — nj)rii = n%S n .fi, since each state can accommodate no more than one 
particle. As the system is in contact with a thermostat, the relative transition rate, 
Tji(nj,ni)/Tij(ni,rij), must also be proportional to a weight factor, e~^ €j ~ €i \ where 
/3 = l/(k B T). 

After setting-up the initial microconfiguration, the simulation proceeds by allowing 
the particles to make random jumps between the single-particle states, with probabilities 
proportional to the transition rates. In what follows we shall denote by P{ ni } the 
probability to find the system in the microconfiguration {rii}. Then the mean occupation 
of a single particle state is 

{".} 
After the system is simulated long enough, the equilibrium is reached and the 
ensemble averages of the thermodynamic quantities do not change in time anymore. 
Once this condition is reached, the equilibrium probabilities, pjj, satisfy the detailed 
balance equation (DBE), 

P{n k } T W k },{n k } = P^n' k } T {n k },W k }- (2) 

for any neighboring states, {n^} = (. . . , rii, . . . , rij, . . .) and {n' k } = (. . . ,n' i7 . . . ,n'j,. . .) = 
(. . . , rii — 1, . . . , rij + 1, . . .). On the other hand, since the system is in contact with a 
heat bath, the equilibrium probabilities should be the canonical probabilities, 

pZ} = e ~^ niei] l Z ^ N ^ ( 3 ) 

where Z(T, N) is the canonical partition function. 

The DBE can be easily checked. If we plug the expression ((3]) into ((2]), with 
T{n' },{n k } = Tji(rij,ni) oc (1 ± nj)rii, we obtain an identity. Therefore the procedure 
described above produces the canonical probabilities when equilibrium is reached. 

3. Simulation of systems which obey intermediate statistics 

3.1. Fractional exclusion statistics 

The calculation of the transition rates in FES systems is not as straightforward as 
in the Bose and Fermi gases. A FES system is essentially an interacting particle 
system, described by a (countable) set of "quasiparticle" quantum numbers, k i; where 
i = 0,1, . . . jU [71 El [T71 HU [in]. For example ID integrable quantum gases of interacting 
bosons and fermions [I~9| |20j [10] and systems described by the Fermi liquid model 
[2T1 l22l li~5l l23l rm [TO] are typical FES systems. For the ID integrable quantum gases 
the quantum numbers are the asymptotic momenta, ki, whereas for the Fermi liquid 
systems the "quasiparticle" energies, ii, are used as quantum numbers. The population 
of the "quasiparticle" states will be denoted by n; or n(ki). The particles are interacting 
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and therefore the values of the quantum numbers may depend (in principle) on the whole 
set {n.i}j=o,i,..., which defines also the quantum state of the system. 

To show how the dependence of kj on {rij}j=o,i,... leads to FES (see also [IE] for 
more details), let us assume that the allowed values of kj, {k;}j =0) i : ..., form a set of 
points in a ^-dimensional (dD) vector space. To build a statistical description of the 
system, we have to calculate the number of configurations in which the particles can 
occupy different one-particle states at constant total energy E and particle number N. 
For this, we divide the dD space of k vectors into elementary, fixed volumes, Vi, each 
containing Gi states and iVj particles - these are our species and should be made no 
confusion from using the same subscripts, i,j, etc. for counting both, the states, |kj), 
and the species, Vi,Gi,Ni. If the particles in the system are bosons, rij may take any 
positive integer value and the number of microconfigurations we have in the system for 
this division into species is 

w »« G '-*» = n wlfr-~i)r (4a) 

If the particles are fermions, Wj = or 1 and the total number of microconfigurations is 

w°>.*})=n m £i m - («) 

Since in our system the values {kj} change with the particle populations, {rij}, 
while the elementary volumes, Vi, are fixed in the k space, this leads to a dependence 
of the GVs on the set {A^} i=0>li ... [9j Q21 [181 HD| • 

To quantify this statement, let us assume that a small change, 5Ni, of iV, produces 
a linear change, 

5Gj = —aji5Ni, (5) 

of the number of states, Gj, for any i and j (not all the a's have to be different from zero). 
The proportionality constants, a^, are the FES parameters defined in the Introduction. 
Using Eqs. (13. ip and (JSJ) we can calculate the grandcanonical partition function, 

Zf&iiGitNi}) = y 52e?Xt N *<r- t *W p ({G i ,N i }) (6) 

{Ni} 

where we assumed that the total energy of the system is E({Ni}) = ^2 i N i €i. This 
is equivalent to ascribing the same energy-the average energy-to all the particles in a 
species. In Eq. (J^D p stands for B or F, depending on whether we have bosons or 
fermions in the system. Using Eq. (jHJ) we can calculate the equilibrium properties of 
the system by maximizing Z FES with respect to the set {N^} pj [10]. In this way we 
obtain for bosons and fermions the systems of equations 
l + fij 
fii 



JJ(1 + hj)-^ = e^'-d (7a) 



1 - hi 



Hi 



Y[(l - nj)^ = e^ e ^ l) , (76) 
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respectively, where n, = Ni/Gi. One should note that in both Eqs. (I7aj) and (ITfej) a^j = 
corresponds to non-interacting bosonic and fermionic systems, respectively. 

To calculate the transition rates, we have to consider separately the systems of 
bosons and the systems of fermions; we shall be interested in transitions in which a 
particle jumps from one species to another, therefore changing the particle distribution 
over the species. 

Let's consider the jump of one particle from the species i into the species j and 
vice- versa, in a system of bosons. For the jump i — y j, we start from a configuration 
of Ni particles and Gi available states in species i and arrive into a configuration of Nj 
particles and Gj states in the species j. The reversed process, j — y i, is exactly like 
i — y j, but in reverse order, i.e. we start from a configuration of Nj particles and Gj 
states in the species j and arrive in a configuration of N particles and Gi states in the 
species i. 

Let's now focus on the process i —y j and assume that the N particles are 
distributed as n.j , n^ , . . . onto the states with the quantum numbers kj ,kj 1; .... 
Analogously, the Nj — 1 particles that exist in the species j before the jump, are 
distributed as n JO , n^ , . . . particles on the states kj , k 7l , . . .. Keeping the notations 
from Section [2], we have the transition rates Tj rnin (nj rn ,n in ) oc (n Jm + l)rii n . Therefore 
the transition rate of a particle-any particle-from the species % to the state k Jm of species 
j is Tj m i(rij m , N) oc {rij m + l)N. Moreover, the transition of a particle from the species 
i onto any state of the species j can be found by summing-up over all the Gj states of 
the species j, leading to TB,ji(Nj, N) oc (Gj + Nj — l)iVj. 

Further, to obtain the total transition rate from the species i to the species 
j, we have to multiply the above result with the number of microconfigurations, 
W Btij = (Gj + Nj - 2)\/[(Nj - l)\(Gj - 1)!] x (G< + N - l)!/[iV,!(^ - 1)!]. In this 
way we finally obtain, up to a common multiplicative factor and the canonical weight, 

e P(ei-ej) 



■ 



T (N IN) ^ + iV ^- 1 ) ! (^ + ^~ 1 ) ! (8) 

T Bjji (Nj - 1, Ni) - {N ._ mG ._ 1)liN ._ mG ._ l)V W 

The DBE then reads, 

e -Wi<*+Vfi-i)<i\T Btii {Nj - 1, N)e^-^ = e-ftM-Vv+^TB.ijiNi - 1, Nj) (9) 

which, by using Eq. ([8]), it can be immediately checked to be an identity. 

A similar argument can be used for transitions in systems of fermions, with the 
difference that the particles cannot arrive on a quantum state already occupied. This 
leads to a transition rate Tpji(Nj, N) oc (Gj — Nj)N. Multiplying again by the number 
of microconfigurations, we obtain 

T Ftji (Nj - 1, N) = (jVj _ 1)! 4_ iVj)!(A ,_ 1)!( G l -iV J )!- ( 10 ) 

Like in the Bose systems, the transition rates (fTOl) satisfy identically the DBE. 
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3.2. Gentile's statistics 

In GS each single-particle state can be occupied by up to £ particles. This can be 
described in a second quantized theory by introducing the creation and annihilation 
operators a k and a k , respectively, of the form [2] 

/ \ / 1 



1 





V2 



\ 



Ok 



v^ 0/ 








V2 



\ 



\ 



o 



[ID 



J 



k x ~ ' k 

The creation and annihilation operators are assumed to commute for different k's. If 
we denote by |0) the vacuum state, then a many-body state of iV = n ko + n kl + . . . 
particles is |n ko ,n kl , . . .) = (4 ) ?lko ( a k 1 ) nkl • • • I )- 

As the system is connected to a heat bath, the eigenstates of the free GS particle 
system are perturbed by external interactions. The transition rate of a particle from 
the state |kj) to the state \kj), Tji, caused by this interaction should be proportional 



to a k .aki- Using (fill) , we obtain Tji oc (ra kj + l)n ki , like in the case of bosons, with the 
extra condition that Tji = if n k > £. These transition rates obey the DBE, like in 
Section El 

The thermal properties of the GS have been calculated in several papers (see for 
example Ref. [2]). We assume that the single-particle states |kj) are eigenstates of the 
Hamiltonian of the (ideal) system, with the eigenvalues e, (or, equivalently, e k J. Since 
only up to £ particles can occupy at same time such a state, the partition function can 
be written as 



£{GS) 



n 



J=0 



jPfa-ei) 



n 






'12a) 



From Eq. fll2a[) we can calculate the grandcanonical thermodynamic potential, 



QiGS) = -k B T\og(Z^) 
and the total particle number, N = d\og(Z ( - GS ^) /d(/3fx) 

y^t ne n/3(fi-ti) 



~(GS) 

n) 



?n/3(n—ei) 



1 



— Y^/i n i ' wriere 

=(5+i)/3(^-m) _ 1 



126) 



:i2c) 



is the average population of the single-particle level i. From the Eqs. (I12q|) . we can 
calculate all the thermodynamics of the system. 

At this point the differences between FES and GS become clear. We can split 
the GS system also into species, like the FES system. Although in both statistics the 
maximum number of particles per species is the same (provided that the species have 
the same dimensions), the number of microconfigurations is different in general. This 
has been calculated in several papers, e.g. [H El El]. If in the species i of a GS system 
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we have, like in FES, Gi states and N t particles, and we denote by Ri = [Gi/(£ + 1)] 
the greatest integer contained in Gj/(£ + 1), then the number of microconfigurations in 
which the particles can be distributed in the species i is 

W (G AH-YV , Y Gi[Gi + Ni-^ + l)r-l]\ ( 

W GS (G h Ni ) - ^(-1) rKGi _ r)m _^ + 1)r] r (13) 

Expression ( lT3"j) is different, in general, from Eq. ( 14 op . although it reduces to it in special 
cases, like e.g. when £ > Gi in ({TBI and a = in (j4aj). 

These differences are reflected in the thermodynamic properties of GS and FES 
systems, a fact which is particularly obvious in systems with constant density of states. 

A transition rate between different species in GS would be rather difficult to write, 
since it depends on how many states in the target species are fully occupied. Therefore 
for GS systems we shall consider just transitions from one state to another, as discussed 
in the beginning of this section. 

The main thermodynamic quantities of GS have been calculated for example in 
Ref. j2]. In the next section we will recover these results by Monte Carlo simulations. 

4. Numerical implementation 

In order to simulate the time evolution of the system and, in the same time, to validate 
numerically the transition probabilities introduced in Section El we take as model system 
the dD ideal gas, contained in a box of linear dimension L, with the density of one- 
particle states 

y{ ' 2m \2ii) V{d/2) K ' 

We split the energy axis into "short" segments, (e*, e i+ i), i = 0,1, . . ., which define the 
particle species. The segments are chosen in such a way that each species contains the 
same number of single particle states, G = J 1+1 g(e)de, for any i. We shall consider 
only the case a^ = aSij . For the convenience of numerical calculations, we require that 
G/a is an integer, while for GS we take G — 1. 

We set the total number of particles in each FES system to be N and we fix the 
energy scale by setting the Fermi energy in each system. For the numerical calculations, 
we have to fix the highest energy in the system, e max , so that the contribution coming 
from the higher energy levels can be neglected in the thermodynamic results. 

In this paper we shall consider only ID and 2D systems, so we define a scaling 
temperature, T , by pA d (T ) = 1, where p = N/L d is the dD particle density and 
A(T) = h/y/2-n-mkBT is the thermal length. This gives 

N 2/d / h 2 \ 1/2 

T *=-FT A: ,d = l,2. (15) 

Concretely, we shall take €p(a) = p{T = 0; a) = a and e max = 30, which would 
permit us to vary the temperature in a range between zero and 2T , without any 
observable effect on the thermodynamic quantities. We set the energy unit to ep(a = 1). 
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In FES systems we cannot use directly the transition rates given in Eqs. flSJ) 
and ffTUl) . because they involve in general too big numbers which cannot be calculated 
numerically The quantities that are interesting for us are actually the relative transition 
rates, T p ji(Nj — l,Ni)/T Pi ij(Ni,Nj — 1). Moreover, in order to obtain a common 
description for both, bosons and fermions, we define Goi = Gi + aiVj, for bosons, 
and Goi = Gi — (1 — a)N t + 1 for fermions. In these notations the number of 
microconfigurations for the species i becomes [7J |25] 

W{G ^ Nl) ~ m(G K - a Ni - 1)! (16) 

for any type of particles. Therefore from now on we shall omit the subscripts B, F, or 
p in our notations. 

To calculate the ratio T p> ji(Nj — 1, Ni)/T Pt ij(Ni, Nj — 1) we take its logarithm and, 
by noticing from Eqs. (jHJ) and (TTU1) that T p ji(Nj — 1, N) = T P)i j(Ni — 1, Nj), we write 

-TjtiNj-^N] 



log 



[TyiNitNj-l)} 



logiTijiNi - 1, Nj)] - \og[T tJ (N t , Nj - 1)] 
dN { dNj 



l0g ^ [l+(l-o)n,l 1 -°(l-«n,)° c _ lif . f {l ' } 

hi 

where Ui = N/Gi and in obtaining the last line we used the Stirling approximation, 
logiV! ps i\nog(iV/e). Therefore, from now on we shall use 

Tji{Nj, Ni) = fn[l + (1 - ^^^-"(l - an 3 ) a , (18) 

up to the relative canonical factor, e^ ei ~ ej > . One can understand the derived transition 
rates as follows: a new state is proposed by moving one particle from species i to species 
j with a step probability which depends on the instantaneous average occupations of 
the two species, i.e. ~ n,[l + (1 — a)n,,] 1_a (l — afij) a , and it is accepted with the usual 
Metropolis probabilities, min(l, exp(/3(e£— 6j)). In this way the dynamics of FES systems 
can be viewed as the dynamics of classical ideal gas, but with new step probabilities 
which stem from a generalized exclusion principle. This picture is consistent with the 
one described in Ref. [5] for Bose and Fermi systems. From Eq. ( TIT]) or (TTBl) it becomes 
obvious that when equilibrium is reached, all the transition rates become equal and 
therefore the DBE is satisfied, as shown in Section 13.11 

In the numerical calculations we shall take Goi = G = 120 for any % and iV = 12000. 
These values are large enough to obtain the analytical equilibrium distribution with 
high accuracy. Nevertheless, a discussion about the species' size is also given below. 
Using these parameters, the considered systems are proved to be very well described 
analytically by continuous models, so that a direct comparison is possible. 

Typical results for equilibration in two-dimensional systems with the FES 
parameter a taking the values 0, 1/4, 1/2, 3/4, 1 are indicated in Fig. fTJ The systems 
are first prepared in their ground state (T = 0i^) configurations and the dynamics 
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10 




2x10 3x10 

step 



Figure 1. Typical equilibration of two-dimensional FES systems, by raising (main 
plot) or lowering (inset) the temperature. The horizontal lines mark the average 
equilibrium energy of each system, calculated analytically. 




Figure 2. Main plot: Analytical (solid lines) and numerical (symbols) data for a two- 
dimensional system with a number of single particle states in a species G = 120, for 
the five indicated temperatures. Inset: Numerical data for T = O.O6X07 with G = 120 
(solid) and G = 1,2,4, 8 (dashed), where the largest deviation corresponds to G = 1. 
Here the symbols have been dropped for convenience. 



towards equilibrium is depicted at T = To. The inset contains the opposite situation 
i.e. the systems are quenched from T to 0.5T . In this case, the initial configuration 
is a random configuration drawn according to the equilibrium distribution, calculated 
analytically at To. In both cases, i.e. either increasing or decreasing the temperature, 
after a number of ~ 10 6 steps, the energy of systems gets close to and starts oscillating 
around the equilibrium value, as calculated analytically. 
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Figure 3. Comparison between FES and GS for semions (a = l/£ = 1/2), for low 
(T = O.lTo) and high temperatures (T = To). Analytical data is represented by solid 
lines and numerical data by symbols. 
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Figure 4. Average total energy vs. temperature (symbols) for d = 1 (a) and d = 2 (b) 
for FES systems with a — 0, 1/4, 1/2,3/4, 1. Inset: Temperature dependence of the 
chemical potential. Solid lines represent analytical calculations using the corresponding 
continuous model. 



Figure [2] shows the average occupation of a two-dimensional system with a = 1/2, 
which is calculated performing time averages after the system has reached equilibrium. 
The system is first very well thermalized for number of A^ = 10 9 steps and then 
the time average is obtained using an equal number of additional steps. We included 
every 1000 th microconfiguration in the time-average, i.e. we considered a total of 10 6 
micro configurations . 

The data is fitted very well by analytical calculations using the corresponding 
continuous model. From Eq. (17)ofRef. [7], one can see that in the degenerate limit, i.e. 
low temperatures, the different distributions cross at the point (E = /xt=o = 0.495, (n) = 
a/4/5), as it is indicated by dotted lines, the small deviation of 0.005 in //t=o, being 
caused by the finite spacing between the energy levels. In FES systems, if one decreases 
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Figure 5. Heat capacity vs. temperature (symbols) for one-dimensional (a) and two- 
dimensional (b) FES systems, and for two-dimensional GS systems (c). Solid lines 
represent analytical calculations using the corresponding continuous model. 



the number of single particle states in a species, G, the average distributions deviate, as 
it is shown in the inset. For the limiting case, G — 1, it is established that the average 
FES distributions for different temperatures in the degenerate limit cross at (n) = 1. 
Therefore we emphasize the high importance of the species' size in obtaining the correct 
FES distribution. 

In Fig. [3] we analyze comparatively the two statistics, FES and GS. For this we 
plot the equilibrium distributions for semions (a = l/£ = 1/2). Although in the high 
temperature limit the distributions become both identical with the classical Maxwell 
Boltzmann distribution, there are visible differences at low temperatures. One particular 
difference is that GS equilibrium distributions in the low temperature limit cross at 
(n) = £/2, which is not the case for the FES systems. 

Next we perform long runs to obtain thermodynamic quantities with good accuracy. 
As before, we take N = 10 9 steps to reach equilibrium and another Ni = 10 10 
steps to obtain the time averages. The results for the average energy and chemical 
potential as function of temperature are depicted in Figs. HI The chemical potential is 
extracted as a fit parameter for the average distribution obtained numerically, in order 
to match the analytical result for the mean occupation. The characteristics related to 
the dimensionality of the systems are clearly visible: for d = 2 case the data is only 
vertically shifted, which is a consequence of the thermodynamic equivalence of systems 
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of constant density of states and any statistics [HI H2J EH US] • 

The heat capacity is calculated from the fluctuations of the total energy, C v = 
({E 2 ) — (E) 2 ) / {k-QT 2 ) . The data plotted in Fig. [5]confirms again the correctness of our 
approach for the proposed step probabilities for FES and GS systems. As expected, 
for FES systems and d = 2, the data for all values a collapses on the same curve. 
By contrast, the thermodynamic behavior of GS systems is no longer equivalent with 
respect to the £ parameter. The thermodynamic equivalence of the 2D systems suggests 
that FES is a more natural extension of Fermi and Bose statistics. 

5. Conclusions 

The present approach introduces a stochastic model for the non-equilibrium dynamics of 
systems which obey generalized statistics, such as fractional exclusion statistics (FES) 
and Gentile statistics (GS). The model has at its core a time- dependent Markov chain 
in the microconfiguration state space, which generalizes previous results from Ref. [5]. 
In the derivation of the stochastic transition probabilities, a division can be made which 
separates them into acceptance probabilities, corresponding to transition probabilities in 
a classical ideal gas (i.e. Metropolis probabilities), and step probabilities, which account 
for the generalized exclusion principle employed. 

The obtained probabilities are tested extensively by Monte Carlo simulations on 
one- and two-dimensional FES and GS systems. The numerical results reflecting 
several thermodynamic quantities overlap very well with the reference data, calculated 
analytically. We point out here the crucial role played by the size of individual 
species in the FES systems, as well as the fundamental differences which appear in 
the thermodynamics properties of the two-dimensional FES and GS systems. 

The Monte Carlo approach to FES could yield new insights in the time evolution 
of finite quantum many-body systems regarded as FES systems, indicating potential 
differences related to out-of-the-equilibrium phenomena. 
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